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OBJECTIVES  OF  THE  RESEARCH  EFFORT 


The  objectives  are: 

•  To  examine/establish  a  novel  method  for  understanding  the  mechanisms  involved  during  non¬ 

destructive  and  destructive  testing  of  (composite)  materials  better. 

•  To  utilize  the  technique  as  a  tool  for  close  cooperation  between  diverse  materials  evaluation 

groups,  i.e.  to  establish  concurrent  analytical,  computational,  non-destructive  and  destractive 
research  on  composite  materials. 

Both  objectives  have  been  scrutinized,  and  the  results  have  been  very  encouraging.  Detailed 
examination  of  the  “transverse”  test  (details  are  given  below)  through  computer  simulations  has 
shown  clearly  the  importance  and  relevance  of  the  first  objective.  Since  the  simulation  process 
makes  possible  the  “visualization”  of  the  mechanisms  undergoing  on  the  interface  subjected  to  load, 
all  groups  listed  under  the  second  objective  have  interest  in  utilizing  the  technique.  For  example, 
as  detailed  in  the  sequence,  it  is  possible,  using  the  technique,  to  investigate  and  visualize  the 
influence  of  residual  stresses.  Since  this  influence  depends  on  the  processing  temperature  as  well 
as  on  the  geometry  of  the  specimen/structure  examined,  groups/individuals  working  on  the 
processing  aspect  of  the  composite  materials  examined  have  interest  in  utilizing  the  simulation 
technique. 


STATUS  OF  THE  RESEARCH  EFFORT 

Soon  after  this  grant  was  awarded,  it  was  decided  to  concentrate  on  the  so-called  transverse  test 
which  aims  at  the  mechanical  properties  of  the  fiber-matrix  interface  in,  mainly,  MMCs.  The  reason 
for  this  decision  was  the  availability  of  unique  and  very  relevant  nondestructive  and  destructive  tests 
performed  at  the  Materials  Laboratory,  Wright-Patterson  Air  Force  Base  [1].  Based  on  [2],  the  test 
configuration  (the  transverse  test  described  in  detail  in  the  sequence)  is  discretized  into  a  lattice  which 
delineates  the  matrix,  the  fiber,  and  the  interface.  The  simulations  provide  further  understanding  of 
the  irechanisms  involved  during  the  relevant  testing.  In  addition,  through  back-analysis,  quantitative 
values  of  the,  homogenized,  interface  properties  can  be  obtained. 

In  addition  to  the  response  of  the  interface  to  mechanical  loading,  the  influence  of  residual 
stresses  have  been  examined  in  detail.  In  an  earlier  study  [3]  the  effect  of  specimen  geometry  on  the 
residual  stress  distribution  was  exarrnned,  mainly  near  firee  surfaces.  The  simulations  have  shown  that 
the  specimen  geometry  as  well  as  the  processing  temperature  are  key  factors  in  determining  the 
interfece  response.  Especially  for  fracture  initiation,  the  residual  stresses  are  more  often  than  not  the 
most  crucial  governing  factor.  As  an  important  consequence,  definition  of  interface  properties 
independent  of  residual  stresses  is  of  dubious  merit.  The  present  study  suggests  strongly  that  the 
interface  cannot  (should  not)  be  examined  in  the  context  of  local  continuum  theories,  but  rather  as 
an  entity  closely  related  to  the  specimen/structure  geometry  it  is  embedded  in  as  well  as  to  the 
processing  temperature.  The  simulation  results  agree  closely  with  the  analytical  solution  [3]  where 
the  effect  of  residual  stress  near  boundaries  is  addressed.  Further,  the  technique  can  be  applied  to  any 
problem  independent  of  the  existence  of  an  analytical  solution.  Detailed  results  are  presented  in  the 
relevant  publications. 
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BACKGROUND 

There  is  a  recent  effort  to  investigate  the  transverse  properties  of  titanium  matrix  composites  (TMC) 
reinforced  with  continuous  silicon  carbide  (SiC)  fibers,  and  this  involves  underetanding  of  the 
dependence  of  the  fiber-matrix  debonding  on  the  residual  stresses,  on  the  fiber-matrix  bond  strength, 
and  on  the  matrix  properties,  under  transverse  loading  conditions  of  the  composite.  Novel  methods 
for  testing,  nondestructively  as  well  as  destructively,  the  damage  evolution  rnechanisms  of  a 
transversely  loaded  sanple  have  been  developed  [1].  In  these  experiments  samples  with  a  single  fiber 
were  used.  Fig.  1. 


Fiber 


Loading 

Direction 


FIG  1:  Schematic  of  Transverse  Test  Configuration  Showing  Orientation  of  Fiber  in 
a  Sample  and  the  Direction  of  Loading 


FIG  2:  Evolution  of  Damage  During  Transverse  Loading  of  a  Composite  Specimen. 
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The  evolution  of  interface  damage  during  the  transverse  test  is  illustrated  schematically  in  figure  2, 
[1]  and  references  cited  therein.  During  the  initial  stages  of  loading,  the  radial  compressive  residual 
stress  (not  near  boundaries)  at  the  interface  decreases  while  the  interface  remains  intact.  As  the 
external  loading  increases  further  progressive  failure  of  the  chemical  bond  occurs,  starting  from  the 
two  sides  in  the  direction  of  loading  and  progressing  to  the  entire  interfacial  debond. 

Figure  3  shows  [1]  the  stress  strain  response  obtained  from  a  typical  transverse  test  specimen, 
and  figure  4  [1]  shows  the  “Reflectivity  Amplitude”  at  various  stages  of  loading,  as  obtained  by 
relevant  ultrasonic  testing.  The  higher  the  reflectivity,  the  more  debonded  the  interface.  Thus,  figures 


FIG  3:  Stress-Stain  Response  of  Specimen  Subjected  to  Transverse  Loading 

3  and  4  show  clearly  that  interface  failure  does  not  occur  simultaneously  nor  homogeneously.  Rather, 
the  process  of  failure  is  highly  heterogeneous  along  the  length  of  the  fiber,  and  no  specific 
“heterogeneity  pattern”  can  be  identified. 

The  numerical  simulations  of  the  transverse  tests,  cf.  Appendices  agree  both  qualitatively  and 
quantitatively  with  the  above  experimental  observations,  described  briefly  herein.  Some  limitations 
of  the  simulation  approach,  at  the  present  stage  of  development,  are  discussed  in  the  next  section  as 
well  as  in  the  appendices. 


LIMITATIONS  AND  FUTURE  WORK 

As  is  shown  in  both  appendices,  A  and  B,  the  numerical  simulations  are  performed  using  a  so-called 
lattice  analysis.  Undoubtedly,  a  three  dimensional  analysis  is  preferable.  This  is  especially  true  for 
simulation  of  the  transverse  test,  since,  due  to  boundaiy  effects  in  the  vicinity  of  the  fiber  edges,  a 
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FIG  4:  Spatial  Along  the  Fiber  Variation  ofR^ectivity  Amplitude  at 
Various  Stages  of  Loading  Indicated  in  Figure  3 

representative  plane  stress  or  plane  strain  analysis  of  the  problem  cannot  be  pursuedAdentified. 
However,  a  fully  three-dimensional  lattice  analysis  is  not  trivial.  This  is  due  to  the  “unaccepta.ble” 
conputer  time  required  for  3-D  lattice  analyses.  The  problem  can  be  partially  resolved  by  using  either 
ultra  fast  single  processor  computers,  i.e.  large  scale  super  computers,  or,  preferably,  parallel 
(massively)  computing. 

Utilmng  advanced  computational  facilities  and  techniques  was  out  of  the  scope  of  this  project 
However,  since  even  the  2-D  analysis  pursued  shows  promising  results,  it  seems  only  natural  to 
extend  the  present  work  to  fiiUy  3-D  analysis.  This  would  alleviate  some  of  the  problems  encountered, 
cf.  Appendices. 


PUBLICATIONS 

G.  ftantziskonis,  Theodore  E.  Matikas,  Prasanna  Karpur,  S.  Krishnamurthy,  and  Leon  Shaw,  Lattice 
Analysis  to  Assess  Fiber-matrix  Interface  Behavior  under  Various  Experimental 
Configurations,  In:  Computational  Mechanics '  95,  S.N.  Atluri,  G.  Yagawa,  T.A.  Cruse,  Eds., 
Vol  2,  pp  2563-2568. 

The  above  publication  is  listed  in  APPENDIX  A 

G.  Frantziskonis,  T.E.  Matikas,  P.  Karpur,  The  Effects  of  Residual  Stresses  on  the  Interface 
Response  in  MMCs,  Simulation  and  Experiments,  paper  under  preparation,  to  be  submitted 
to  Composite  Structures  or  other  relevant  Journal. 

The  major  results  for  the  above  publication  appear  in  APPENDIX  B 
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INTERACTIONS  AND  ADVISORY  FUNCTIONS 

The  first  of  the  two  papers  listed  above  was  presented  by  George  Frantziskonis  at: 

Conputational  Mechanics '  95,  International  Conference  on  Computational  Engineering  Science,  July 
30, 1995,  Hawaii. 

The  project  was  completed  under  continuous  interaction  with  the  Materials  Laboratory,  Wright 
Patterson  Air  Force  Base,  Ohio.  The  PI,  George  Frantziskonis,  interacted  continuously  with  scientists 
there,  cf.  following  names.  He  visited  WPAFB  in  July  1995  to  continue  the  cooperation  and 
interacted  with: 

Dr.  Thomas  Moran,  WL/MLLP 

Dr.  Theodore  Matikas,  UDRI,  on  site  at  WL/MLLP 

Dr.  Prasanna  Karpur,  UDRI,  on  site  at  WL/MLLP 

One  graduate  Research  Assistant,  Mr.  Gang  Hong  has  been  supported  from  this  project.  Part  of  the 
results  presented  herein  are  to  appear  in  his  Dissertation,  scheduled  for  final  examination  during  Fall 
1996. 
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1.  INTRODUCTION 

This  work  utilizes  a  novel  procedure  for  obtaining 
mechanical  properties  of  the  fiber-matrix  fiber 'bn^^ag^matrix  yield  and/or 

simulates  actual  expenmenis  in  ‘**‘“*’  1  oroceduie  was  implemented  for  the 

cracking,  and  interlace  fmlure.  In  a  J'  the^fraamentation  test  for  metal  matrix 

following  commonly  *  h™u”?nd'ouilout  tests^for  MMCs  as  well  as  ceramic 

composites  (MMCs):  b)  'h®  “e  test  configuration  is  discretized  into  a  lattice 

matnx  composites  In  *®  and  the  can  be  found  in  (1].  The 

which  delineates  the  matnx.  the  fiber.  mechanisms  involved  during  the 

simulations  provide  funher  '^“derstandmg  ot  the  rnechan^  J 

retevaa.  teanng.  In.  SrpS  wt  fira.  dascrib. 

snnnlau/na  for  d.  so<^ad  nansvc« 

nf  rhe  same  maicrial  composition. 


2.  FRAGMENTATION  TESTS 

The  .ensue  load  Is  applied  in  *n  direcuon  of  the 

i£ioSSei»7pps=s 

.0  ac&ve  avenge  ftagmenudon 

be.w«.n  nber  -d  .minx  is 
modeled  IS  ai  Shase.' Two  valoes  am  necessary  .0  descr.be  .m  elasne  bn.de- 
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response  -  Ore  in.en.hase  sum«s  c«ffic«n.  S  (modulns  'f  “ 

Stress  coefficient  F  matnx.  SiC  (SCS-6)  fiber  composite  were 

fragmentation  the  values  for  Ti  fragmentauon  pattern 

idcniificd  as:  S=860  MP^fm.  F-  example,  the  fragmemauon  patiem 

SS-Sr^rabT^'vlal  m™'Jsens,.ive  (w.Upn  .00%  changes,  .0  .he  v^ue  o. 
S,  but  depends  strongly  on  F. 


2.  PUSHOUT  AND  PULLOUT  TESTS 


Such  tests  are  commonly  fie^  wUh  Us  endsT^^^ 

composites.  The  specimens  cont^n  “  ^mg^*  J  ^  compression  unul  it  slides  out  of  the 
the  pushout  test  an  indenter  loads  the  tiber  comp  ^ 

matnx.  The  pullout  test  ‘Vro^o^^  'he  following:  (a)  smoil. 

fiber.  For  metal  matnx  compos  external  load  application  side,  for  both  the 

stable  cracks  fomi  at  the  mentioned  above,  tension  and 

pushout  and  pullout  tesa.  Ih)  |j.j,iate  at  the  side  opposite  to  the  external  load 
shear  dominated  interface  ^  that  initiated  at  the  bottom,  funher  external 

applicauon,  (0  after  the  amst the  size  of  the  "hole  ’  at  the  bottom 
load  is  required  for  the  ^ho‘e  mtertac^  ,he  effects  of 

S  S"S“£“p«nnnu  Dc«.Ucnn  b.  found  in  d»  d). 

3.  TRANSVERSE  TESTS 

The  .nns«r«  .«.  which  is  used  ” 

usuaUy  performed  using  single  ^o^el  co  p  /  hemicalbonding.  which 

suited  for  interface  the  fiber  mUtenals  dunng 

depends  on  the  chemical  h  samole  as  well  as  in  a  multi-ply  composite 

or^essing.  remains  the  same  in  a  single  fiber  s^pie  as  wen  ciculate  in  the 

snJ  (b,  .h.  ^htoaJ  sd«s«  in  dK 

Lgle  fiber  sample,  thereby  making  it  fet^ible  to  account  lo 

msSoUng  of  U..  exporimenud  «ul.s 

the  fiber-matnx  interface.  In  ^e  to  lowing  modelmg  studies, 

from  recent  transverse  tests  ^cn^e  re  “  ^cre  processed  by  the 

Monofilament  composite  samples  fwthesyx^^me  y 

foil-fiber-foil  technique  wherein  two  Ti-6A1-4V  sh«B  wim  a  smgi^ 

SCS-0)  between  them  were  hot  orocessing  condition  clearly 

hours.  The  interfacial  gl^Mrbon  coadng  and  the  matru 

?^.T.n"?.rdSSnTi;c"rSr^^^|;u— 

ro“d“  g^-Sne'^Srs^r^^^^  with  I  fiber  axis  perpendicular  to  the  loading  axis  of 
the  samples  as  shown  in  Figure  la. 

3.1  Experimental  Approach  o„t  using  a  micro-straining  stage  built  in  the 

WUMaST==.™.T31.Thn.W^^ 

“r’S“siSct'^un 

using 


1  Ar  ^ach  stcD  of  lo^dinc,  the  fibcr-msirix 

samples  was  done  in  janjple  under  dial  load)  in  a  pulse-echo  mode 

interface  was  imaged  ^Sident  on  the  surface  of  the  samples  at  an 

using  a  focused  ultrasonic  beam  which  was  mci 

angle  of  24®.  an  angle  polarized  shear  waves  were  incident  on  the 

material.  Hence.  to  the  transducer  and  the 

fiber-matnx  ‘"“^^ice.  The  shear  waves  to  characterize  the  interface 

reflection  coefficwnt  from  the  ^  modeling  of  the 

fracture.  The  modeled  based  on  theoretical 

reflection  coefficient  " “d  ampli^^^^ 

predicuons  and  expcnmental  Jn  o  P  ^^^rate  the  maximum  reflected 

4V/SCS-6  single  fiber  oompiwite  samp  ^  [g^gls  of  the  sample 

corresponding  ultrasonic  shear  a  labeled  ‘A’  in  Figure  lb  corresponds  to 

labeled  ‘A  through  K’  in  the  Figum  TS^^^cmenVof  the  mechanical  loadinV  of  the 
the  fiber-matrix  interface  f.^t  few  points  of  interface  that  fracture  at 

sample.  Reflectivuy  .mage  ®  Lno^SeSra^^^  to  red  in  the  color  scale) 

about  350  MPa  (shown  with  the  miwinmm  ^pmuae  ^ 

These  points  of  interface  f^mte  are  .j.  the 

places  at  the  center  of  the  fiber.  The  reflect  y  ag^^^  .^s 

progression  of  interfacial  damage  as  important  conclusion  drawn 

eVre  interface  has  been  f«ctujed  at  ab^t 

from  using  in-situ  SBR  <„  jp  fracture  almost  instantaneously  once  a 

assumption  that  the  entire  interface  is  “*“‘y  ,  .  existence  of  a  weak  diffusion 

sufficiently  high  stress  level  reached  be«u^^^^ 

bond  (61  (as  contrasted  to  pSnts  at  a  low  load,  to  the 

the  debonding  progresses  from  a  nu  3^0  MPa  as  shown  in  Figure 

entire  interface  over  a  finite  range  of  pp  redistribution  of  stress  along  the  interface 

w‘h\c?'ocS’rrdue‘^  mJhe”ir?pagadon 

3.2  Lattice  Modeling  of  Transverse  Tests  three-dimensional.  Thus  a  simulation 

The  test  configurations  menuoned  atove  ^  Umiiations 

should  reflect  the  3-D  effects.  ^^l^g-jensively  in  [1].  For  the  fragmentation  and 
from  such  a  "simplificaiion  annlied  oarailel  to  the  fiber,  the  limitations  from  a 

pullout,  pushout  tests,  since  the  ^  although  a  3-D  simulation  process 

2-D  simulations  may  not  be  of  8r«at  “npo^^  •  2-D^imulations  have  been  the 
would  be  preferable.  The  reason  P*^°  .  ®  d  i  e  thermally  induced  residual 

following:  (a)  there  are  several  issues  ^ 

stresses,  the  impheauons  of  “5^iI“nnVin^  tte  ^  fiberl^d  matrix,  etc.,  before  a 

on  interface  properties.  d  (hi  for  3-<iiinensionai  analysis,  only  problems 

fully  3-D  simulation  process  is  employed,  (b)  for  3-^mcnsio^  ^  y  ^ 

with  rather  coarse  discretization  can  be  important 

‘h&».  c«J= 


^  TO 


...  r  .r-vp'^/  ',  '  V''^ 


progressively,  (d)  the  matrix  starts  yielding  only  after  "complete"  failure  of  the  interface, 
after  that  stage  the  matrix  starts  yielding  and  develops  a  "shear  band"  type  failure  as 
shown  in  Figure  4,  (e)  the  external  load  increase  needed  to  advance  interface  failure  from 
initial  to  its  final  failure  state  is  not  negligible. 


Figure  1  (a)  Experimental  Configuration  Showing  Transverse  Orientation  of  the  Fiber  in 
a  Sample  and  the  Direction  of  Loading,  (b)  In*Situ  Ultrasonic  SBR  Imaging  of  an 
Embedded  SCS-6  Fiber  in  Ti-6A1-4V  Matrix  During  Various  Stages  of  Transverse 
Loading. 


Figure  2:  (a)  Two-dimensional  configuration  of  the  transverse  test.  The  bars  in  the  lattice 
assigned  interface  propenies  are  shown.  The  fiber  is  inside  the  "ring,"  and  the  rest  of  the 
configuration  is  matrix,  (b)  Interface  crack  pattern  from  two-dimensional  simulation  of 
the  transverse  test.  The  corresponding  external  load  level  is  shown  by  an  arrow  in  Fig.  3. 


OOFY  ATAIUBUI 


TO  DTIC  DOES  NOT  PERMIT  FULLY  LEGIBLE  REPRODUCTION 


In  our  simulations  S*d?f&uU  to°defiM  anextemal  "stress"  that 

progressed  up  to  a^ut  520  ^  ‘  nn^ental  configuration.  Despite  this  difficulty 

corresponds  to  interface  failure  fortes  ex^n^  ^  considered  a 

,n  order  to  examine  the  response  along  the  <hrecuo^  Performed  simulations  (2-D).  as 
"failure  stress"  of  420  MPa  I^^SleTto^^e  ^  and  the  y-axis  coincides 

shown  in  Figure  5a.  where  now  Ae  ,  interface  elements  are  subjected  to 

with  the  external  load  axis,  k  these,  s  that 

ite  same  ares,  to  w.  chose  .0  „,f„„ider.  »U1  be  impenam  for 

residualsresses.  winchrhe  ptont  wy  0  inienace  crek 

mterface  fracture  Here  tatetto  fracture  .s 

tom  utmrface  tour,  .nttiauou  to  its  complete  tocture  ui 


Figure  3:  The  external  suess  Vs  external  straia  respous.  as  calculated  from  m.  two 
dimensional  simulation. 


Ftgure  4t  Simulatiou  reults.  Failure  (y»‘'^) 
interface  faUure.  The  corresponding  external  strain  is  0.02. 

Thus,  for  the  !he  dhecdon^of^ke  is  rather 

the  fiber  axis  is  non-brittle,  while  cr^k  8  .  different  directions.  In  other  words, 

brittle.  This  shows  the  strong  interactions  ^ong  te  different^c 

three-dimensional  effects  arc  important  .  configuration  of  Figure  2.  It 

values  can  be  obtained,  however,  by  2-0  y  .gH  iicinz  interface  values  as 

n“e  again  not«l  to  dm  simulattons  She"  suits  of  the  2-D 

?Sii:Ss  comp^J'wto  (averaged  over  .be  f.ber  dlrectton, 


COP^ 


riu-^ 


txperimcnal  respoi«.  Averaging  “  'tocnon.  We'S^cS?^^^ 

cimiiiations. 


ngure  5. 

«  ma«x. »,  m^r/ace  cracic  paa«n  »™ 

after  cm^  iniciaiion. 
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Two  major  improvements  of  the  work  presented  in  Appendix  A  have  been  investigated; 

-  Detailed  study,  using  the  lattice  based  simulation  process,  to  determine  the  residual  stress 

distribution,  i.e.  before  any  external  mechanical  load  is  applied  on  the  specimen. 

-  Examination  of  the  fiber-matrix  interface  properties  including  the  residual  stress  distribution.  In 

other  words  this  second  “improvement”  calls  for  identification  of  the  apparent  interface 
properties  which  include  the  effect  of  residual  stress. 

Before  presenting  and  discussing  the  results,  the  mechanical  properties  used  for  the  titamum  alloy 
matrix  and  for  the  SiC  fiber  are  given  in  the  following  table.  The  processing  temperature  was 

considered  to  be  900®C. 


Density 

gr/cm^ 

Young’s 

Modulus 

GPa 

Shear 

Modulus 

GPa 

Yield 
(matrix)  or 
Failure 
(fiber) 
MPa 

Poisson’s 

Ratio 

Thermal 

Expansion 

Coeff 

lOYC 

Ti-6A1-4V 

45 

0.32 

SiC  SES-6 

179 

3500 

0.15 

1  3.99 

As  mentioned  above  the  analysis  was  restricted  to  two-dimensions.  Thus,  first  the  response  of  the 
specimen  far  from  the  fiber  edges,  i.e.  for  an  infinite  sized  structure  was  analyzed  by  using  a 


FIG  Bl:  The  First  2-D  Configuration  Analyzed,  (a)  80x80  Lattice  where  the  Bars 
Assigned  Fiber  Properties  are  Shown,  (b)  the  Interface  Bars  are  Shown. 


configuration  such  as  that  of  Rg.  Bl.  The  lattice  sizes  varied  firom  80x80  to  100x100.  The  dimensions 
of  each  specimen  are  0.0015x0.0015m,  the  fiber  diameter  is  150pm,  and  the  interface  properties  are 
considered  (normalized/divided)  per  unit  of  its  thickness,  as  explained  in  a  lucid  fashion  in  [2].  Thus, 
in  the  sequence  all  interface  properties  are  given  as  values  per  thickness.  In  the  configuration  shown 
in  Fig.  Bl,  the  interface  thickness  is  0.0015/80  m,  that  of  one  lattice  spacing. 

The  matrix  is  considwed  elastic-plastic,  the  elastic  part  being  linear  and  the  hardening  being 
also  linear  (elastic,  linear  hardening).  In  other  words,  the  matrix  is  allowed  to  yield  and  develop  high 
strains.  The  fiber  is  considered  linearly  elastic  and  brittle.  After  the  failure  stress  is  reached  the  load 
(i.e.  the  load  of  a  firactured  bar  in  the  lattice  analysis)  is  reduced  to  zero.  The  interface  is  also 
considered  to  be  linearly  elastic  and  brittle.  Then,  as  discussed  in  detail  in  [2]  the  problem  reduces  in 
determining  the  elastic  and  failure  properties  of  the  interface.  This  is  done  by  back  analysis  of  the 
simulations  described  below  and/or  the  values  determined  from  back  analysis  of  tests  other  than  the 
transverse  are  used.  For  example  for  Ti-6A1-4V  matrix,  SCS-6  fiber  composite  the  interface 
properties  were  determined  in  [2]  by  back  analysis  of  the  fragmentation,  the  pullout  and  the  pushout 
tests. 

Figure  B2  shows  six  (6)  views  of  the  residual  stresses  before  any  external  load  is  applied,  i.e. 
due  to  the  processing  terrperature  difference.  As  can  be  seen,  the  residual  stresses  are  quite  high,  i.e. 
as  compared  with  the  yield/failure  stresses  for  the  matrix/fiber.  Note  that  this  configuration.  Fig.  Bl 
corresponds  to  an  iii^tely  long  fiber,  thus  boundary  effects  are  excluded.  Thus  the  interface  is 
mostly  subjected  to  compressive  stress.  The  elasticity  properties  used  for  the  interface  in  obtaining 
the  stress  distribution  shown  in  Fig.  B2  are  those  obtained  by  back  analysis  of  the  firagmentation, 
pullout  and  pushout  tests  [2].  The  so-called  modulus  coefficient  (Young’s  modulus  divided  by 
interface  thickness)  is  857  MPa/pm.  In  [2]  the  interface  failure  strength  coefficient  was  determined 
to  be  11  MPa/pra  In  the  transverse  test  the  thickness  of  the  sample  is  0.004m,  thus  the  fact  that  the 
behavior  near  the  edges  is  actually  different  (boundary  effects)  is  important.  Thus,  failure  is  governed 
by  the  interplay  of  boundary  and  bulk  (homogeneous)  interface  response.  It  was  determined  that  the 
value  of  11  MPa/pm  for  the  interface  failure  strength  coefficient  needed  to  be  adjusted  for  the 
(apparent)  inclusion  of  boundary  effects  to  13.65  MPa/pm.  Once  again,  problems  of  such 
“adjustments”  would  not  appear  if  a  fully  three  dimensional  analysis  was  performed.  Figures  B2a  and 
B2b  show  a  3-D  plot  of  the  residual  stresses,  for  the  configuration  of  Fig.  Bl.  They  show  that 
a  low  tensile  stress  develops  in  the  matrix,  while  a  “high”  con^ressive  value  develops  at  the  interface 
and  in  the  fiber.  Hgures,  B3c-f  show  contour  plots  of  the  same  stress  distribution,  Oqq,  included  herein 
for  the  sake  of  clarity. 

The  residual  stresses  influence  the  process  of  failure  of  specimens  significantly.  For  the 
configuration  of  Fig.  Bl  partial  cracks  of  small  extend  develop  in  the  interface  during  the  early 
straining  stages.  However,  these  cracks  did  not  propagate  and  cover  only  a  small  fraction  of  the  total 
interface.  Before  further  interface  cracking  occurs,  a  large  part  of  the  matrix  yields.  Figure  B3  shows 
a  typical  interface  fracture  “pattern”  at  some  stage  of  loading.  From  that  point  on,  the  interface  failure 
is  rather  rapid  and  the  total  interface  fails  at  an  external  strain  of  about  3%. 

In  order  to  study  the  distribution  and  the  effects  of  residual  stresses  near  the  fiber  ends,  the 
configuration  of  Rg.  Bl  is  inadequate.  Again,  a  three  dimensional  analysis  would  be  preferred.  Here 
we  “bypass”  a  fully  3-D  analysis  by  studying  the  2-D  configuration  like  that  of  Fig.  B4.  Figure  B5 
shows  different  views  of  the  residual  stress  o^y.  The  high  tensile  stresses  near  the  edges  are  evident 
This  agrees  with  experimental  observations  which  call  for  tensile  microcracks  being  present  at  the 
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interface,  near  the  edges  of  the  specimen.  Thus,  the  interface  is  already  (before  the  external  straining 
is  applied)  cracked  near  the  specimen  edges.  This  inidal  cracking,  together  with  the  residual  stresses 
in  the  remaining  of  the  interface  influence  the  way  the  interface  cracks  under  external  straining  of  the 
sample.  As  shown  in  Fig.  B6  subsequentiy  to  the  initial  cracking,  interface  cracks  develop  near  the 
two  sarrqjle  edges  as  well  as  in  the  center.  This  agrees  with  the  experimental  observations  described 
above. 

Note  that  the  configuration  of  Fig  B4  is  not  representative  of  the  “true”  response,  due  to  the 
limitations  of  a  2-D  analysis.  Despite  this,  note  that  the  stress-strain  response  obtained  from  the  lattice 
analysis  agrees  very  well  (quantitatively)  with  the  experimental  one,  that  of  Fig.  3,  page  4. 
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FIGB4:  The  second  2-D  configuration  Analyzed,  (a)  40x100  Lattice  where  the  Bars  Assigned 
Fiber  Properties  are  Shown,  (b)  the  Interface  Bars  are  Shown. 
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1.  INTRODUCTION 

This  work  utilizes  a  novel  procedure  for  obtaining  quantitative  information  on  the 
mechanical  properties  of  the  fiber-matrix  interface  in  composite  materials.  The  method 
simulates  actual  experiments  in  detail,  including  fiber  breakage,  matrix  yield  and/or 
cracking,  and  interface  failure.  In  a  recent  study,  the  procedure  was  implemented  for  the 
following  commonly  performed  experiments:  (a)  the  fragmentation  test  for  metal  matrix 
composites  (MMCs);  (b)  the  pushout  and  pullout  tests  for  MMCs  as  well  as  ceramic 
matrix  composites.  In  the  simulations,  the  test  configuration  is  discretized  into  a  lattice 
which  delineates  the  matrix,  the  fiber,  and  the  interface.  Details  can  be  found  in  [1].  The 
simulations  provide  further  understanding  of  the  mechanisms  involved  during  the 
relevant  testing.  In  addition,  through  back-analysis,  quantitative  values  of  the, 
homogenized,  interface  properties  can  be  obtained.  In  this  paper,  we  first  describe 
simulations  of  the  pushout/pullout  and  fragmentation  tests  for  a  titanium  matrix,  silicon 
carbide  (SiC)  fiber  composite.  Relevant  interface  properties  are  evaluated  by  simulating 
the  former  test.  Using  these  values,  we  study  the  response  of  both  test  configurations, 
and  then  compare  the  numerical  results  with  actual  experimental  data.  Further,  we 
combine  recent  experimental  results  with  relevant  simulations  for  the  so-called  transverse 
test  of  the  same  material  composition. 


2.  FRAGMENTATION  TESTS 

The  tensile  load  is  applied  in  the  direction  of  the  fiber.  Load  is  transferred  from 
the  matrix  to  the  fiber,  which  at  some  point  of  loading  breaks.  Further  loading  results  in 
the  fiber  breaking  into  successively  smaller  fragments  until  the  fragments  become  too 
short  to  enable  further  increase  in  the  fiber  stress  level.  Since  it  is  the  fiber-matrix 
interface  that  actually  delivers  the  load  from  the  matrix  to  the  fiber,  the  fragmentation 
"pattern”  depends  strongly  on  the  properties  of  the  interface.  It  is  not,  however,  the 
interface  that  solely  governs  the  fragmentadon  pattern.  Relevant  simuiadons  [1]  suggest, 
in  addition  to  the  importance  of  the  interface  properties  discussed  in  the  next  paragraph, 
that:  (a)  the  fiber  fragmentadon  pattern  depends  strongly  on  the  volume  of  the  matrix 
present,  (b)  the  matrix  hardening  modulus  Emh  influences  the  fragmen^don  pattern,  (c) 
the  overall  load-displacement  response  of  the  specimen  is  pracdcally  insensitive  to  the 
interface  properties,  (d)  it  is  pracdcally  impossible  to  achieve  average  fragmentadon 
lengths  of  the  order  of  or  lower  than  the  fiber  diameter. 

In  the  numerical  analysis,  the  transition  region  between  fiber  and  matrix  is 
modeled  as  an  "interphase."  Two  values  are  necessary  to  describe  its  clastic  brittle- 
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„spo.s.  -  *e  =f?S 

stress  coefficient  F  ^ 'if  L%v  matrix.  SiC  (SCS-6)  fiber  composite  were 

fragmentation  “st-  the  values  fwTi-6AMJ^  fragmentation  pattern 

identified  as:  S=860  MP^fm  F-ll  MP^im.  u  is  fragmentauon  pattern 

SSne^wd^^  abo^^^^^^^^^  ?l]"stse°nsitive  (within  100%  changes)  to  the  value  of 

S,  but  depends  strongly  on  F. 


2.  PUSHOUT  AND  PULLOUT  TESTS 


such  tests  are  commonly 

composites.  The  specirnens  contain  “  ^  j  compression  until  it  slides  out  ot  the 
the  pushout  test  an  indenter  lof  s  he  tensile  load  is  applied  on  the 

matrix.  The  pullout  test  is  s‘™lar.  with  the  diHejen^  following:  (a)  small, 

fiber.  For  metal  matrix  Pf  “f  thrextemal  load  application  side,  for  both  the 
stable  cracks  form  at  the  interface  near  ’^he  e  mentioned  above,  tension  and 

p„shou,  and  pullom  »s.s,  (b)  •(«'  *  =  .he  ««mnl  load 

shear  dominated  interface  L  thal  initiated  atV  bottom,  futthet  external 

application,  (c)  after  the  amst  of  J f  ^^e  size  of  the  ''hole  "  at  the  bottom 

load  is  required  for  the  ‘f ‘f^,^f„\'°iefonnatk,n  response,  and  (f)  the  effects  of 
Sffi  SSSm'be  itnpodant.  Details  can  be  fonnd  in  the  litetatnte  (1|. 


3.  TRANSVERSE  TESTS 


The  transveme  test  tehtch  is  n»d  “ 

usually  performed  using  single  f^her  mo  P  chemical  bonding,  which 

suited  for  interface  evaluation  because,  (a)  materials  during 

depends  on  the  gjg  fjfjej-  sample  as  well  as  in  a  multi-ply  composite 

processing,  remains  the  same  in  a  smg  •  ,  are  relatively  easy  to  calculate  in  the 
panel,  and  (b)  the  residual  stresses  residual  stresses  in  the 

Lgle  fiber  sample,  thereby  m^ng  exclusions  about  the  stress  at  fracture  of 

modeling  of  the  test  ro  hen^e  m  describe  some  experimental  results 

the  fiber-matnx  interface,  J?  J°  results  from  the  lattice  modeling  studies, 
from  recent  transverse  tests  *  experiments  were  processed  by  the 

Monofilament  S  wX  single  SiC  fiber  (SCS-6  or 

foil-fiber-foil  technique  wherein  “f  960  OC  at  17  MPa  for  1.5 

SCS-0)  between  them  were  hot  a  processing  condition  clearly 

hours.  The  interfacial  carbon  coating  and  the  matrix 

rnrdt'tan^SpSSScna  «ilb  *c  fiber  axis  perpendichirh  .o  the  lo.dihg  axis  of 
the  samples  as  shown  in  Figure  la. 

3. 1  Experimental  Approach  ucinB  a  micro-straining  stage  built  in  the 

Transveree  tensile  tests  ^re  c^  g  3^  j^at  the  ultrasonic 

WLAlatenals  Direc™^^^^  at  different  stress  levels.  An 

S  uKonic  "“""bve 
using 


1  in  mrrPTn^ntal  steos  At  each  step  of  loading,  the  fiber-matrix 

anele  of  24°  an  angle  between  the  first  and  the  second  cntical  angles  of  the  m 

s;ro„c„®.=^ 

J'jJrS  ^Sr^^vS-sSo"!!  ™  J  >0  c>lita»  .he  maxim™  mnec.eh 

Signal  from  the  fiber  to  correspond  to  the  test  of  a  T1-6A1- 

4wscs?™sirrrrni;s’?a;;prs=“ 

S'iSngimsi.upR^ 

=1 

ih'i  This  ranee  is  dependent  mainly  on  the  redistribution  of  stress  alo  g 

^d^he  Image  of  the  fiber  significantly  reduced  as  the  open  interface  cracks  closed. 

m  “Igh'SSSS  above  am  tee^onal.  Thus  a  stadata 
u  M  rhf*  n  effects  In  fll  the  simulations  were  two-dimensional,  the  limitations 

S“stphto^  in  [1].  For  the  fragmentation  and 

pullout,  pushout  tests,  since  the  l^oad  is  applied 

2  D  simulations  mav  not  be  of  great  importance,  although  ^  3  D  P 

would  be  orefLSilI  The  reason  for  performing  2-D  simulations  have  been  the 
following-  ^(a)  there  are  several  issues  to  be  resolved,  i.e.  therm^ly  induced  residu^ 
stresses,  frie  iinplications  of  assuming  a  homogenized  'nte/phase  effect  of  tempemtur^ 
on  interface  properties,  imperfections  in  the  interface,  fiber,  and  naatnx,  betore  a 
fuilv  3-D  simulation  process  is  employed,  (b)  for  3-dimensional  analysis,  o  y  p 
S  iSer  "dlLtization  caL  solved  within  reasonable 

the  transverse  test,  the  limitations  of  a  2-dimensional  analysis  may  be  importanh 
Sow^ver.  certain  issues  can  be  examined  by  analyzing  a  configuration  like  the  one 

shown  2b  shows  the  crack  pattern  obtained  from  the  u 

level  shown  by  an  arrow  in  Figure  3.  The  following  have  been  idenufied.  (a)  before  *e 
Sx  Slds.  L  interface  starS  to  fail,  (b)  interface  faUure  is  not  symmetnc  w«th  ^ 
to  either  the  loading  axis,  or  the  axis  perpendicular  to  the  loadir^  direction,  (c) 

Interface  failure  is  locally  brittle.^he  final  stage  of  interface  failure  is  reached 
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progressively,  (d)  the  matrix  starts  yielding  only  after  "complete"  failure  of  the  interface, 
after  that  stage  the  matrix  stans  yielding  and  develops  a  "shear  band"  type  failure  as 
shown  in  Figure  4,  (e)  the  external  load  increase  needed  to  advance  interface  failure  from 
initial  to  its  final  failure  state  is  not  negligible. 


Figure  I  (a)  Experimental  Configuration  Showing  Transverse  Orientation  of  the  Fiber  in 
a  Sample  and  the  Direction  of  Loading,  (b)  In-Situ  Ultrasonic  SBR  Imaging  of  an 
Embedded  SCS-6  Fiber  in  Ti-6A1-4V  Matrix  During  Various  Stages  of  Transverse 
Loading. 


Figure  2:  (a)  Two-dimensional  configuration  of  the  transverse  test.  The  bars  in  the  lattice 
assigned  interface  properties  are  shown.  The  fiber  is  inside  the  "ring,"  and  the  rest  of  the 
configuration  is  matrix,  (b)  Interface  crack  pattern  from  two-dimensional  simulation  of 
the  transverse  test.  The  corresponding  external  load  level  is  shown  by  an  arrow  in  Fig.  3. 


In  our  simulations  interface  failure  initiated  at  about  339  MPa  extern^  sttess  ^ 
progressed  up  to  about  520  MPa.  Thus  it  is  difficult  to  define  an  external  str«s  that 
LtrLponds  to  interface  failure  for  this  experimental  confiprauon.  Despite  Jfficulty 
in  order  to  examine  the  response  along  the  direcuon  of  the  fiber,  we  have  considere 
"failure  stress"  of  420  MPa  and  with  this  value  we  performed  simulauons  (2-D),  as 
sSwn  .to  now  ■]»  x-axin  in  pamllei  to  thn  flben  to  to  y-™.nc.to 

with  the  external  load  axis.  In  these  simulauons.  all  interface  eleiMnts  ^  s^jected  to 
the  same  stress,  thus  we  chose  to  initiate  fracture  at  the  right  ^ 

residual  stresses,  which  the  present  analysis  does  not  consider,  will  be  important  for 
interface  fracture  along  the  fiber  direction.  Figure  ob  shows  the  interlace  crack 
propagation  at  some  point  of  external  load  applicauon.  Here  interface  fract^  is  bnttie 
(the  external  stress  increase  from  interface  fracture  imuauon  to  its  complete  fracture  is  not 
appreciable). 


Figure  3:  The  external  stress  Vs  external  strain  response  as  calculated  from  the  two- 
dimensional  simulation. 


Figure  4:  Simulation  results.  Failure  (yielding)  of  the  ductile  matrix,  occurring  after 
interface  failure.  The  corresponding  external  strain  is  0.02. 

Thus,  for  the  Unnsverse  test,  interface  (overall)  failure  on  the  plane  transverse  to 
the  fiber  axis  is  non-brittle,  while  crack  propagation  in  the  direction  of  the  fiber  is  rather 
brittle.  This  shows  the  sttong  interactions  among  the  different  directions.  In  other  words, 
three-dimensional  effects  are  imponant  for  the  transverse  test.  Mean  interface  failure 
values  can  be  obtained,  however,  by  2-D  analysis  using  the  configurauon  of  Figure  2.  It 
is  once  again  noted  that  the  simulations  were  performed  using  interface  values  ^ 
determined  from  the  back  analysis  of  the  fragmentation  test.  Thus  the  results  of  the  2-D 
simulations  can  be  compared  with  the  mean  (averaged  over  the  fiber  direction) 


experimental  response.  Averaging  is  necessary  since  the  3-D  interactions  may  render 
interface  failure  highly  non-uniform  along  the  fiber  direction.  We  are  currently 
process  of  correlating  experimental  and  simulation  re'sults,  and  on  performing  actual  j-D 
simulations. 


Figure  5:  (a)  Along  the  fiber  axis  two-dimensional  configuration  of  the  t^sverse  test. 
The  bars  in  the  lattice  assigned  interface  properties  are  shown.  The  fiber  is  inside  the  two 
parallel  lines  and  the  rest  of  the  configuration  is  matrix,  (b)  Interface  crack  pattern  soon 
after  crack  initiation. 
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